A sparse matrix formulation of model-based ensemble Kalman filter
نویسندگان
چکیده
Abstract We introduce a computationally efficient variant of the model-based ensemble Kalman filter (EnKF). propose two changes to original formulation. First, we phrase setup in terms precision matrices instead covariance matrices, and new prior for matrix which ensures it be sparse. Second, split state vector into several blocks formulate an approximate updating procedure each these blocks. study simulation example computational speedup approximation error resulting from using proposed approach. The is substantial high dimensional vectors, allowing run on much larger problems than can done with In introduced block negligible compared Monte Carlo variability inherent both procedures.
منابع مشابه
The Ensemble Kalman Filter : theoretical formulation and practical implementation
The purpose of this paper is to provide a comprehensive presentation and interpretation of the Ensemble Kalman Filter (EnKF) and its numerical implementation. The EnKF has a large user group, and numerous publications have discussed applications and theoretical aspects of it. This paper reviews the important results from these studies and also presents new ideas and alternative interpretations ...
متن کاملA deterministic formulation of the ensemble Kalman filter: an alternative to ensemble square root filters
The use of perturbed observations in the traditional ensemble Kalman filter (EnKF) results in a suboptimal filter behaviour, particularly for small ensembles. In this work, we propose a simple modification to the traditional EnKF that results in matching the analysed error covariance given by Kalman filter in cases when the correction is small; without perturbed observations. The proposed filte...
متن کاملA Matrix-Free Posterior Ensemble Kalman Filter Implementation Based on a Modified Cholesky Decomposition
In this paper, a matrix-free posterior ensemble Kalman filter implementation based on a modified Cholesky decomposition is proposed. The method works as follows: the precision matrix of the background error distribution is estimated based on a modified Cholesky decomposition. The resulting estimator can be expressed in terms of Cholesky factors which can be updated based on a series of rank-one...
متن کاملA Local Ensemble Kalman Filter for the Ncep Gfs Model
This paper outlines the basic concept and mathematical formulation of the Local Ensemble Kalman Filter (LEKF, Ott et al. 2002 and 2003) data assimilation scheme. Some important properties of the scheme are illustrated by numerical experiments with the Lorenz-96 model. Initial implementations of the LEKF on the National Centers for Environmental Prediction Global Forecast System (NCEP GFS) are d...
متن کاملResampling the ensemble Kalman filter
Ensemble Kalman filters (EnKF) based on a small ensemble tend to provide collapse of the ensemble over time. It is shown that this collapse is caused by positive coupling of the ensemble members due to use of one common estimate of the Kalman gain for the update of all ensemble members at each time step. This coupling can be avoided by resampling the Kalman gain from its sampling distribution i...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
ژورنال
عنوان ژورنال: Statistics and Computing
سال: 2023
ISSN: ['0960-3174', '1573-1375']
DOI: https://doi.org/10.1007/s11222-023-10228-0